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Abstract 

The Configuration Interaction (CI) method using orbitals centered on the nucleus has recently 
been applied to calculate the interactions of positrons interacting with atoms. Computational 
investigations of the convergence properties of binding energy, phase shift and annihilation rate 
with respect to the maximum angular momentum of the orbital basis for the e + Cu and PsH bound 
states, and the e + -H scattering system were completed. The annihilation rates converge very slowly 
with angular momentum, and moreover the convergence with radial basis dimension appears to be 
slower for high angular momentum. A number of methods of completing the partial wave sum are 
compared, an approach based on a AX j = a(J + f)- n + &(J+f)-( n+1 ) form (with n = 4 for phase 
shift (or energy) and n = 2 for the annihilation rate) seems to be preferred on considerations of 
utility and underlying physical justification. 
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I. INTRODUCTION 



In the last few years there have been a number of calculations of positron-binding to atoms 
IE IE 0, IH IE IE IE IE EE, U| and positron-scattering from atoms and ions [IE EE 3 
using orthodox CI type methods with large basis sets of single center orbitals. Besides the 
calculations on atoms, there have some attempts to calculate the properties of positrons 
bound to molecules with modified versions of standard quantum chemistry methods 2> 3 
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One feature common to all these CI type calculations is the slow convergence of the 
binding energy, and the even slower convergence of the annihilation rate. The attractive 
electron-positron interaction leads to the formation of a Ps cluster (i.e. something akin to 
a positronium atom) in the outer valence region of the atom |E IE llil Tgj ] . The accurate 



representation of a Ps cluster using only single particle orbitals centered on the nucleus 
requires the inclusion of orbitals with much higher angular momenta than a roughly equiva- 
lent electron-only calculation 0,IE HE 21 1- F° r example, the largest CI calculations on PsH 



and the group II positronic atoms have typically involved single particle basis sets with 8 
radial functions per angular momenta, £, and inclusion of angular momenta up to £ max = 10 
[E IE IE]- Even with such large basis sets, between 5-60% of the binding energy and some 
30-80% of the annihilation rate were obtained by extrapolating from the £ max = 10 to the 
^max = oo limit. 

Since our initial CI calculations on group II and IIB atoms |5j, |y, |7j , advances in com- 
puter hardware mean larger dimension CI calculations are possible. In addition, program 
improvements have removed the chief memory bottleneck that previously constrained the 
size of the calculation. As a result, it is now appropriate to revisit these earlier calculation 
to obtain improved estimates of the positron binding energies and other expectation values. 
However, as the calculations are increased in size, it has become apparent that the issue of 
slow convergence of the physical observables with the angular momenta of the basis orbitals 
is the central technical issue in any calculation. 

Whilst it is desirable to minimize the amount of mechanical detail in any discussion 
(so as not to distract from the physics), the ability to draw reliable conclusions from any 
calculation depends crucially on the treatment of the higher partial waves. For example, 
the treatment of the higher partial waves in separate calculations by Saito Hj has already 
been shown to be flaw ed |22j while the present work exposes the defects in the methods of 
Gribakin and Ludlow [14]. The present work, therefore, is solely devoted to an in-depth 
examination of the convergence properties of mixed electron-positron calculations. 

In our previous works, [E IE S EE EE| , a relatively simple solution to this problem was 
adopted. In effect, it was assumed that the successive increments to any observable scaled 
as an inverse power series in £ max . This approach does have limitations as do the approaches 
adopted by other groups ^EEI> 3 IE|- I n the present article, we examine the convergence 



properties of positron binding calculations upon PsH and e + Cu, and a positron scattering 
calculation upon the e + -H system with respect to angular momentum and the dimension 
of the radial basis sets. The limitations of existing calculations are exhibited, and some 
improved prescriptions for estimating the variational limit are introduced and tested. 
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II. EXISTING METHODS OF PERFORMING THE ANGULAR MOMENTUM 
EXTRAPOLATION 



A. The nature of the problem 

The positron-atom wave function is written as a linear combination of states created by 
multiplying atomic states to single particle positron states with the usual Clebsch-Gordan 
coupling coefficients ; 

\*;LS) = ^2a d (LM^mjlLM^iSiMs^jlSMs) 

x $i(Atom; LiSi^jiTo) . (1) 

In the case of a single electron system, e.g. H, Q^Atom; LiSA is just a single electron 
wave function, i.e. an orbital. For a di-valent system, <&i(Atom; LiSi) is an antisymmetric 
product of two single electron orbitals coupled to have good Lj and Si quantum numbers. 
The function <pj{ro) is a single positron orbital. The single particle orbitals are written as a 
product of a radial function and a spherical harmonic: 

0(r) = P(r)Y tm (r) . (2) 

The radial wave functions are a linear combination of Slater Type Orbitals (STO) 23] 
and Laguerre Type Orbitals (LTOs). Most of the time the radial functions are LTOs, 
the exceptions occurring for single electron states with angular momenta equal to those 
of any occupied core orbitals. Since the Hartree-Fock core orbitals are written as a single 
combination of STOs, some of the active electron basis is written as linear combinations of 
STOs before the switch to a LTO basis is made. The LTO basis 0, 0| has the property that 
the basis can be expanded toward completeness without having any linear independence 
problems. 

The present discussion is specific to positronic systems with a total orbital angular mo- 
mentum of zero. It is straight-forward to generalize the discussion to states with L > 0, but 
this just adds additional algebraic complexities without altering any of the general conclu- 
sions. 

For a one electron system, the basis can be characterized by the index J, the maximum 
orbital angular momentum of any single electron or single positron orbital included in the 
expansion of the wave function. 

For two electron systems, the L = configurations are generated by letting the two 
electrons and positron populate the single particle orbitals subject to two selection rules, 

max{e ,e 1 ,e 2 ) < J , (3) 

mm{£ 1 ,£ 2 ) < . (4) 

In these rules l§ is the positron orbital angular momentum, while i\ and l 2 are the angular 
momenta of the electrons. The maximum angular momentum of any electron or positron 
orbital included in the CI expansion is J. The other parameter, L int is used to eliminate 
configurations involving the simultaneous excitation of both electrons into high i states. 
Double excitations of the two electrons into excited orbitals are important for taking electron- 
electron correlations into account, but the electron-electron correlations converge a lot more 
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quickly with L mt than electron-positron correlations do with J. Calculations of the positronic 
bound states of the group II atoms and PsH |]| |(| showed that the annihilation rate changed 
by less than 1 % when L- m % was varied from 1 to 3. The present set of calculations upon PsH 
were performed with L int = 4. Further details about the methods used to perform the 
calculations can be found elsewhere 0, 0| . 

Various expectation values are computed to provide information about the structure of 
these systems. All observable quantities can be defined symbolically as 

w J = E AJi > ( 5 ) 

L=0 

where AX J is the increment to the observable that occurs when the maximum orbital 
angular momentum is increased from J— 1 to J, e.g. 

AX J = (X) J - (Xy- 1 . (6) 

Hence, one can write formally 

oo 

(x)°° = (x) J + J2 AxL ■ ( 7 ) 

L=J+1 

The first term on the right hand side will be determined by explicit computation while the 
second term must be estimated. The problem confronting all single center calculations is that 
most expectation values, (X) J converges relatively slowly with J and so the contribution 
of the second term can be significant. A sensible working strategy is to make J as large as 
possible while simultaneously trying to use the best possible approximation to mop-up the 
rest of the partial wave sum. 



B. Existing extrapolation techniques and their limitations 

One of the first groups to confront this issue and attempt a solution was the York Uni- 
versity group of McEachran and Stauffer. They per formed a series of polarized orbital 
calculations of positron scattering from rare gases The decrease in energy 

when the target atom relaxed in the field of a fixed positron was used to determine the 
polarization potential as a function of the distance from the nucleus. The polarized orbital 
method implicitly includes the influence of virtual Ps formation (within an adiabatic approx- 
imation), and this means that slow convergence can be expected. McEachran and Stauffer 
found that the scattering observables, namely the phase shift, and the Z e s annihilation pa- 
rameter, converged slowly with J, the largest angular momentum of the polarized electron 
orbital set used to represent the adjustment of the atomic charge cloud in the field of the 
positron. They took this into consideration by assuming their polarization potential scaled 
as J~ p and the polarized orbital scaled as J~ q at large J. They found p ~ 3.8 and q ~ 1.8 
for the rare gases at J ~ 12. 

The recent Cl-type calculations of Mitroy and collaborators also used an inverse power 
relation of J~ p , to complete the partial wave sum 0, [l2|. In this case, the observables, 
e J , 5 J , r J and Z^ s were extrapolated. This contrasts with the polarized orbital calculations 
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where the polarization potential and polarized orbital were extrapolated. The value of p 
was given by 

/AX J ' 1 \ / fJ-1 



p=ln U^J/ ln l— J • (8) 

while the constant factor is 

A x = AX J J P . (9) 

The correction factor was then evaluated by doing the sum Y1T=j+i Ax/ L p explicitly with 
an upper limit in the thousands. 

Gribakin and Ludlow [Hf applied perturbation theory and the ideas of Schwartz [2£| |3(j 
to determine the asymptotic behavior of the partial wave increments to the binding energies, 
phase shifts and annihilation rates of positron-atom systems. This work is lar gely derived 
from previous work on the partial wave expansion of two electron atoms |3lL l32l I33L l33 |. 
They determined that the binding energy E, annihilation rate T, phase shift 5, and collisional 
annihilation parameter Z e R obey 

AE J = (E) J — (E) J ^ 1 ~ (10) 
at j = ( T y _ ^y-i _ Bt (H) 



(J 



A5 J = (5) J - (5) J -i ~ — %- (12) 



AZ J %S = (Z eS ) J - (Z*)'- 1 ~ jj^- 2 (13) 
These expressions are merely the leading order terms of a series of the form 

To perform the actual extrapolation during a calculation of positron-hydrogen scattering, 
Gribakin and Ludlow did a fit to calculated values at J = 9 and J = 10 with the 
formulae 

00 

5 = (5r = (sy+ £ -—^ 

L=J+1 1 ^ ^ 



J cfT 



^e ff >°° = (^eff) J + £ 

L=J+1 ^ + 2 



«<^> J + ^T' ( 16 ) 

and so determined Z e g and 5. They used the approximate identities in eqs. (fTHjl and (fT7)j) 
rather than explicitly evaluating the infinite sum. The identities appear to have been derived 
as an approximation to the §^ +l (L + ±)~ 2 dL integral. However, in equating the sum to 
the integral they implicitly assume a rectangle rule representation of the integral which is 
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in error of 5-10% for J G [7, 10] (the net effect of this is that Gribakin and Ludlow state 
that the increments decrease as B/(L + |) n but actually assume a B/L n decrease when 
evaluating the J — > oo correction). A better approximation to the series is obtained by 
using a mid-point rule to represent the integral. Doing this leads to 



« . (17) 

/ /+i (L + i)" (n-l)(J+l)-i 1 ) 

This approximation is accurate to 0.1% for n = 2 and J = 7. 

It will be shown that a more serious problem with the Gribakin and Ludlow methodology 
is that eqs. (JT5j) and (|16|) cannot reveal whether the calculated are deviating from the 
expected asymptotic form. For example, successive increments to either the phase shift or 
Z e Q have usually decreased more slowly with J (for J ranging between 10 and 18) than 
indicated by eqs. (fTUj) -(fT3 |) 0, 0, [l^, 131] • Instead of having p = 2 or p = 4 the successive 



increments often gave slightly smaller values for p. The approach adopted by Gribakin and 
Ludlow is insensitive to these deviations. 

Saito has investigated the structure of the PsH and the Ps-halogen systems with the CI 
method 0, [nj. A Natural Orbital (NO) truncation algorithm based on the energy was 
used to reduce the dimensionality of the secular equations, thus making calculations on 
the heavier halogen atoms viable. Besides using the inverse power series, Saito used the 
functional form 

AX J = io-^sio^+t , (18) 

to complete the partial wave sum for the annihilation rate. This function was not based 
on any physical principles, and its usage was justified on the grounds that the increments 
were decreasing faster than J~ 2 . However, it has been suggested that the annihilation rate 
increments were decreasing too quickly because the dimension of the radial basis used in 
the Ps-halogen calculations was simply too small j22|. So the rationale behind the usage of 
eq. (fTHj) is questionable. 

Some mention must be made of the difficulties associated with the slower convergence 
of the annihilation rate. Consider the PsH system, a calculation with J = 9 gave 72% of 
the total annihilation rate |5j. If one doubled the size of J, then eq. (fTTJl suggests that the 
explicit calculation would only recover 86% of the total annihilation rate. And it would take 
a calculation with J « 250 to recover 99% of the annihilation rate. The situation is even 
more sobering when one considers that the annihilation rate converges faster for PsH (since 
it is the most compact) than for any other positron binding system. 



III. COMPARISON OF EXISTING AND NEW APPROACHES TO THE PAR- 
TIAL WAVE EXTRAPOLATION 

A. The different alternatives 

In order to expose the strengths and deficiencies of existing approaches, very large calcu- 
lations have been performed on three mixed electron-positron systems. These are the e + -H 
scattering system for the t = partial wave, and the bound PsH and e + Cu systems. It 
will be seen that the typical calculations on these real-world systems do not agree perfectly 
with the leading order asymptotic form given by Gribakin and Ludlow, i.e. eqs. (|10|) - (|13j) . 
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Accordingly, six different extrapolation methods for determining the J — > oo correction were 
tested. These were: 

Method p. The successive increments to all quantities are assumed to obey an 

X J = , (19) 

law with the exponent p determined from eq. (120(1 . The value of p derived from three 
successive calculation of X J ~ 2 ,X J ~ 1 and X J is given by 

(Note, in previous works we have used a J~ p series [H H, H, EH). The notations pe, pr, 
ps and pz are used to denote the exponents derived from the partial wave expansions of 
the energy, annihilation rate, phase shift and Z e g The discrete sum over L in (Q is done 
explicitly up to J = 200. The remainder of the sum is then estimated using eq. (|17|) . 

Method p av . This is based on Method p. Three successive calculations for (J— 2), (J— 1) 
and J are once again used with eqs. (fT§j) and (J2Uj) to determine an initial estimate of p . 
Then, p is set to the average of po and the expected value of either 2 or 4. This method 
makes an admittedly crude attempt to correct method p in those cases where ps and pz 
were significantly different from 4 and 2 J3j. Once p has been fixed, eq. (|T9"j) can then be 
used to determine A x and the discrete sum over L in ((7j) is done explicitly up to J = 200. 
The remainder of the sum is then estimated using eq. (J17)) . 

Method GL. The relations eq. (J15)) and (J16)) are assumed to be exact. The two largest 
values of (X) J are used to determine (Z c g)°° and (S)°°. This method mimics the procedure 
adopted by Gribakin and Ludlow 

Method I. The functional form 

AX J = — (21) 
(J+|)« 

is assumed to apply and the AX J increment are used to determine Bx- The exponent n 
is set to 2 for the annihilation rate and 4 for the energy or phase shift. The discrete sum 
over L in ((7|) is done explicitly up to J = 200 and beyond that point eq. (|T7|) is used. This 
method has similarities with the GL method. 
Method II. The functional form 

AA ' J = (J+4)» + (Tt|f+ t ' (22) 

is assumed to apply and the AX J and AX J_1 increments are used to determine Bx and 
Cx- The second term in eq. ()22j) comes from 3rd-order perturbation theory (3ll . 32, |33j]. 



The exponent n is set to 2 for the annihilation rate and 4 for the energy or phase shift. The 
discrete sum over L in (J7J) is done explicitly up to J = 200 and beyond that point eq. (|T7jl 
is used. 

Method III. The functional form 

a X J = Bx i Cx i Dx 

(J + !)« + (J+i)-+ 1 + (J + ^+ 2 ' 
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is assumed to apply and the AX J , AX J_1 and AX J ~ 2 increments are used to determine 
Bx, Cx and Dx- Other particulars are the same as those for Methods I and II. 

Method S. The functional form adopted by Saito, eq. (jTHj) is used. The parameters a, 
(3 and 7 are determined from AX J ~ 2 , AX J ~ 1 and AX J . Then the series is completed by 
summing to J = 2000. 



B. The e + -H scattering system 

The CI-Kohn method has already been used to generate phase shift and annihilation 
parameter data for the e + -H scattering system New calculations with a radial basis 



set (i.e. the number of LTOs per £) of increased dimensionality have been done in order to 
minimize the influence of the radial basis upon any conclusions that are drawn. The present 
investigation examines s-wave e + H scattering at k — 0.4 a^ 1 . 

The largest calculation for e + -H included a minimum of 30 LTOs per I with additional 
LTOs included at small t. Special attention was given to the I = 1 positron basis since 
this channel is responsible for the long-range polarization potential. The dimension of the 
LTO basis was 80 in this case. All radial integrations were taken to 729 on a composite 
Gaussian grid. The earlier calculations of Bromley and Mitroy [12] with a minimum of 17 
LTOs per I will be presented for comparison. Table H] gives the e + -H phase shift and Z e g for 
s-wave scattering at k — 0.4 1 up to J = 12. The larger basis will be referred to as basis 
2 while the older basis will be named basis 1. 

First of all, the values of p$ and pz derived from the earlier 0| and present CI-Kohn 
calculations are shown together in Figure ^ as a function of J. This figure tests whether 
eqs. (|T2"j) and (fE^) describe the behavior of a real-world calculation. 

Neither ps or pz are within 1% of the expected asymptotic value at J = 12 and both are 
approaching the expected asymptotic limit from below. The two calculations give almost 



TABLE I: Results of CI-Kohn calculations for s-wave e + -H scattering at k = 0.4 <2q . The 
variational and close-coupling data in the last three rows are taken from calculations which have 
basis functions which explicitly depend on the electron-positron distance. The J — > oo limits were 
taken at J = 12 using the various extrapolation methods as described in the text. The value of 
p were determined at J = 12 using eq. 1)20(1 . The previous CI-Kohn estimates ^3] include the 
J —* oo correction as determined at that time. 



J 


(5} J (radians) 




(z cS y 





-0.1992097664 




0.4529537597 


1 


-0.01002215705 




1.018467362 


2 


0.05339800596 




1.466777767 


3 


0.08175604707 




1.798430751 


4 


0.09626735469 




2.043817341 


5 


0.1043829191 




2.228559914 


6 


0.1092338304 




2.370688616 


7 


0.1122907225 




2.482386967 


8 


0.1143023638 




2.571897852 


9 


0.1156750002 




2.644890677 


10 


0.1166408472 




2.705331457 


11 


0.1173386269 




2.756059864 


12 


0.1178543969 




2.799146079 




J — > oo limits 






P 


3.6248 




1.9583 


Method p 


0.1200652 




3.34020 


Method p av 


0.1199020 




3.32823 


Method GL 


0.1196691 




3.29464 


Method I 


0.1197593 




3.31675 


Method II 


0.1199484 




3.32750 


Method III 


0.1199456 




3.30092 


Method S 


0.1198834 




3.21989 




Other calculations 






CI-Kohn J -» oo [12] 


0.1198 




3.232 


Optical potential [3Ji, j£6j 


0.1201 




3.327 


Variational a [14, 37, 38] 


0.1198 




3.407 


Close Coupling [39, 40] 


0.1191 




3.332 


a The variational results of van Reeth et. al. [ij Ell are only given in tabular form in 








exactly the same ps while the larger calculation tends to give smaller values of pz with the 
difference becoming greater as J increases. The increasing difference between pz for the two 
calculations suggests that a converged calculation of needs an increasingly larger radial 
basis as J increases. This point is addressed in more detail later. In most of the calculations 
we have performed, the values of v derived from eq. (|2(jp have been slightly smaller than the 
expected value at J ~ 10 0, 0, @, E3|- We also note that in all the calculations we have 
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FIG. 2: The extrapolated J — > oo limit of the phase shift (in radians) 5 for e + -H scattering at 
k = 0.4 Qjq clS cl function of J. The horizontal solid line shows the phase shift of Bhatia et.al. 35]. 
The phase shift without any J — > oo correction is given by the curve labeled (5} J . 



so far performed, the values of p increase steadily (once the broad features of the physical 
system have been achieved) as J increases. While the asymptotic increments to 5 and Z e g 
do not agree exactly with eqs. (fE?|) and (JT3J), the observed trends do appear to be consistent 
with their derived limits. 

Figures El and H3 show the behavior of the extrapolated S and Z e g as a function of J for 
some of the different extrapolations. Table|l]gives estimates of S and Z e s using the calculated 
values at the largest possible J value (i.e. 12) to determine the J — > oo corrections. 

There is one problem in interpreting the results of Figures 121 and H3 and Table |U The 
exact value of Z e g is imprecise at the level of 2-3%. The old calculation of Bhatia et.al. 
using Tij co-ordinates gives 3.327 [H, |36|, the close coupling calculation by Ryzhikh and 
Mitroy using the T-matrix method give s 3.332 00, while the variational calculation of 
van Reeth et.al. gives 3.407 H 1371 . l38j . There is also some scatter in the estimates of the 



phase shifts, but the degree of difference between the Bhatia et.al. and van Reeth et.al. 
phase shifts is only 0.3%. (the T-matrix phase shift is only expected to have an accuracy of 
about 1% 0). 

Figure 121 shows that inclusion of the J — > oo correction leads to greatly improved estimates 
of (5)°°. In terms of their impact, the methods belong to 3 operational classes. Firstly, 
Method p consistently gives the largest values of (S)°°. Fixing ps to a set value at a finite 
J inevitably results in the J — > oo correction being overestimated. For example, using 
Ps = 3.343 at J = 8 to fix p$ for all J results in increments to A5 J that do not decrease 
quickly enough. Methods I and GL, on the other hand, tend to underestimate the size of the 
J — > oo correction since p$ = 4 is fixed prior to the sequence of AS J increments achieving 
the expected (J + i) -4 form. 

Those methods which attempt to allow for deviations from the leading order behavior, 
namely Methods II, III and p av approach the expected J — > oo limit much earlier. Indeed, 
their estimates of (S)°° differ by less than 0.5% at J = 6. Table [I] reveals that these 3 
methods give (S)°° estimates that differ by less than 0.1% at J = 12. Of the 3 alternatives, 
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FIG. 3: The extrapolated J — > oo limit using the different methods to complete the partial wave 
series for the e + -H at k = 0.4 clq 1 a function of J. The horizontal solid line shows the Z P fi of 
Bhatia et.al. |3^ | while the horizontal dotted line shows the Z e g- of van Reeth et.al. 



the 3-term asymptotic series, namely Method III, seems to possess the best convergence 
properties. Method S also appears to give a reasonable estimate of {5)°° when J > 6. 

The tabulated estimates of (5) °° at J = 12 reflect the discussions in the above paragraph. 
Method p gives the largest phase shift while Methods GL and I give the smallest phase 
shifts. The maximum difference between any of the phase shifts is only 0.3% since the net 
effect of the J — > oo contribution to the phase shift is only 2%. 

Figure El for Z e g shows some features in common with Figure 121 Once again, application 
of a J — > oo correction is seen to give much improved estimates of the (Z e $)°° limit. Method 
p also gives the largest estimate of Z e g. Method S gives values of (Z eff )°° that are generally 
the smallest, and Table reveals that it gives a value that is 0.1 smaller than any of the other 
approaches at J = 12. It is not surprising that a method based on a fitting function with no 
physical justification performs so poorly, and its reasonable estimate of J — > oo correction to 
the phase shift can be regarded as a numerical coincidence. No further discussion of Method 
S will be made although values are reported in the later tables for reasons of completeness. 

A more detailed analysis and discussion of Figure H3 cannot be made until the convergence 
properties of the underlying radial basis are exposed. 



C. Convergence properties of the radial basis 

In addition to converging very slowly with J, the annihilation rate also converges slowly 
with respect to the number of radial basis functions since the actual wave function has a 
cusp at the electron-positron coalescence point. A previous CI investigation on helium in 
an £ = model indicates that the electron-electron 5-function converged as 0(iV~ 5 / 2 ) where 



N is number of Laguerre orbitals [41| . And it has also been demonstrated that the relative 
accuracy of the electron-electron 5-function increment for a given size radial basis decreased 
as L increased |42^ . 

Some sample ratios can be used to illustrate these points. The ratio R J Z compares the 
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1.08 
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FIG. 4: The ratio of the increments to (5) J and (Z e g) J (refer to ea. (|24j) ) for e + -H scattering at 
k = 0.4 an 1 as a function of J for the 17 and 30 LTO calculations. 



two calculations of AZ^ for e + -H scattering by determining the ratio for basis 1(17 LTOs) 
and basis 2(30 LTOs). It is defined as 



(AZ, 



R ' z = v "~ff° . (24) 



(AZ, 



effA7 



A similar ratio, Rj can be defined for the increment to the phase shifts. A plot of these 
two ratios is given in Figure 01 while Table |H] lists values of Rj and R J Z for some selected J 
values. Figure |U clearly demonstrates the AZ e J ff converges more slowly that the phase shift 
increments at large J. First, the annihilation rate is more sensitive to the size of the radial 
basis than is the phase shift. Second, the higher partial waves are more sensitive to the size 
of the radial basis than the lower partial waves. For example there was a 4.0% increase in 
AZ| ff between the basis 1 and basis 2 calculations while there was a larger 7.2% increase in 
AZgg. However, the increase in A<5 12 was only 1.7%. 

Since the lack of completeness in the radial basis has the largest impact at high J, it will 
obviously affect the J — > oo correction. For example, Method II gives AZ^ 12 = 0.4695 for 
basis 1. For basis 2, the correction of AZ^ 12 = 0.5284 is significantly larger. 

The implications of these results can be seen by consideration of the Method II plot of Z^ 
depicted in Figure EJ This achieves a maximum value of 3.41 at J = 5, and then decreases 
until it is 3.328 at J = 12. The question to be addressed is whether the decrease from J > 6 
is due to convergence properties of the calculation with respect to J or the convergence 
properties of the basis with respect to N, the number of orbitals per £7 Although there are 
seven estimates of Z e s that lie between 3.30 and 3.34 in Table HI we believe that the true 
value of Z e ff lies closer to 3.407 (the value of van Reeth et.al.) than to 3.327 (the value of 
Bhatia et.al.). 

This interpretation is supported by a crude estimate of the variational limit deduced from 
Figure EJ The assumption is made that the AZ^ increments converge as 0(iV~ 5//2 ) |4l| . 
Consequently, one deduces that the plotted R J Z ratios comprise some 57% of the necessary 
correction to the variational limit. The variational limit for any increment then estimated to 
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TABLE II: Ratios of the partial wave increments to the energy (or phase shift), annihilation rate 
(or Z e ff) taken from the two different basis sets used for each system. 
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be AZ c J ff x (1 + 0.43 x (R J Z — l)/0.57). Applying this correction to the data in Table |U gives 
(Z e g)°° = 3.410 when Method II is used to estimate the J — > oo limit (the actual correction 
of AZ^ff 12 of 0.578 was about 9% larger than the basis 2 value). Determination of the 
variational limit for the phase shift can also be done by assuming that the A5 J increments 
converge as 0{N- 7 / 2 ) |4l|. In this case, the plotted Rj ratios comprise some 76% of the 
correction to the variational limit and thus the final estimates of the phase shift increments 
would be A5 J x (1 + 0.24 x (Rj - l)/0.76). The Method II phase shift only increased by 
1.3 x 10~ 5 radian giving (5) 00 = 0.12008. 

One point from Table warrants special attention. The 3-term asymptotic series, namely 
Method 3, gives a smaller (Z e s)°° than Method I! This seems ridiculous given that p = 1.9583 
at J = 12. The tendency for the AZ e J ff increments to be systematically underestimated 
results in the corruption of the Bz, Cz and Dz coefficients extracted from the 3-term fit 
and renders Method III unreliable for determination of (Z e ff)°°. 

Even though data for only one system has been presented, it is possible to make some 
general comments about the performance of the difference methods since analysis of the e + Cu 
and PsH data will confirm these conclusions (they are also compatible with the results of 
large basis CI calculations of He [42| ) . 

Since the ps and pz exponents tend to be smaller than 4 and 2 respectively at finite J, 
Method p has an inherent tendency to overestimate the J — ► oo corrections. Application of 
this approach in the past has not resulted in any gross errors since the problems associated 
with fixing p at values less than 2 and 4 tend to cancel out the errors associated with a 
radial basis of finite size p |7J. This method should only be applied in situations when the 
asymptotic form of the expectation value under investigation is unknown. 

Method I generally underestimates the J — > oo correction. It gives a useful estimate of 
the J — > oo correction and should mainly be applied to give rough estimates for low precision 
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FIG. 5: The exponents (pe and pr) for 2 different CI calculations of e + Cu as a function of J. 



calculations. Method GL can be regarded as a variety of Method I that happens to give 
inferior J —>■ oo corrections. 

Methods II and p av were seen to give J — > oo corrections that were close to each other 
once the calculation reached a certain value of J. Method II should be be preferred since it 
is founded in correct asymptotics. 

Method III seems to give the earliest reliable estimate of the phase shift. However, it 
should not be applied to the annihilation rate unless the radial basis is substantially larger 
than the present basis. Method III should only be applied in situations where the underlying 
partial wave increments have an accuracy of better than 1% and in addition the increments 
should vary smoothly and not exhibit fluctuations. 



D. The e + Cu ground state 

Table HTT1 gives the e + Cu binding energy and annihilation rate as a function of J up to 
J = 18 for the calculation with 25 LTOs. The table also includes values from a calculation 
with the fixed core stochastic variation method (FCSVM) [1,|43||. The FCSVM basis includes 
the electron-positron coordinate explicitly and is very close to convergence. The FCSVM 
calculation uses a slightly different model potential so it is not expected that the CI energy 
and T should be exactly the same. Figure displays pe and pr versus J for two different 
CI calculations of the e + Cu ground state. One plot is derived from the earlier calculation 
of Bromley and Mitroy 0| which included a minimum of 15 LTOs per £- value (basis 1). 
The present calculation (basis 2) is much larger with a minimum of 25 LTOs per £ value 
(note, more than 25 LTOs were included for £ = 0, 1 and 2 since these make the largest 
contribution to the energy and annihilation rate). 

The plots of pe and pr against J for e + Cu are similar to the plots of ps and pz for 
e~-H scattering. Both exponents are generally smaller than the expected asymptotic limits 
but steadily increase as J increases. The actual value of pr at J = 18, namely 2.030, was 
marginally larger than the expected asymptotic limit of pr = 2. The estimates of (e)°° and 
(r)°° as a function of J are shown in Figures El and 
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FIG. 6: The e + Cu binding energy as a function of J. The different curves use different algorithms 
to estimate the J — * oo correction as discussed in the text. 

Table UTTl reveals an 8% increase in T when compared with the earlier CI calculation value 
of T, namely 0.474 xlO 9 s _1 |8[- This is a consequence of the bigger radial basis used in the 
present work. The selected values of listed in table |H] reveal an 8.7% increase in Ar 16 for 
basis 2 compared to basis 1. It is expected that further increases in the radial basis would 
eventually lead to a pr that was smaller than 2.0 at all J = 18. Again it is noticed that 
(and R J E ) increase with increasing J. 

The estimates of the annihilation rate in table IIHI are all very close together (with the 
exception of Method III). This occurs because pr = 2.030 is very close to the expected value 
of 2.0. The variations between the different approaches are largely concerned with taking 
care of the deviations from the pr = 2 behavior, and with a minimal deviation at J = 18, 
one should expect minimal differences between the final results. Method III is the least 
accurate (discounting Method S) since it is the most susceptible to the inaccuracies in Ar J 

Figure [7| shows that Method p systematically overestimates (r)°° at smaller values of 
J. Method I, on the other hand generally gives the smallest estimates of the (T)°° and is 
consistently too small at lower J. 

Method III obtains a reasonable estimate of (e)°° the quickest. Beyond J > 12 the 
Method III binding energy does decrease slightly. This may be due to slower convergence 
of the radial basis at higher J. Methods II and p av give roughly equal estimates of (e) 00 for 
J > 10, and the spread between the Methods II, III and p av estimates of (e)°° is only 0.2% 
at J = 18. 

Method p gives the largest estimate of (ft) 00 for all J shown in Figure 13 while Methods I 
and GL give the smallest values (with Method GL once again being worse that Method I). 
It is worth mentioning that Methods II, III and p av are all roughly constant after J = 16 
while Methods I, GL and p are still increasing or decreasing. 

In summary, the totality of information in Table 11111 and Figures El El an d is very 
reminiscent of the situation for e + -H scattering and is consistent with the conclusions derived 
from the analysis of e + -H scattering. 
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FIG. 7: The e + Cu annihilation rate (in units of 10 9 s 1 ) as a function of J. The different curves 
use different algorithms to estimate the J — > oo correction as discussed in the text. 

E. The PsH ground state 

The best two electron system for validation purposes is the positronium-hydride (PsH) 
system since its properties are very well known as a result of previous investigations [44| 
Ea . l46l I47I ]. The stochastic variational method (SVM) expectation values listed in Table 
IIVI are taken from a new calculation with 1800 ECGs. The energy of this wave function, 
E = —0.7891674 Hartree is the best PsH energy so far reported [44]. 

The orbital basis used in the present calculation was about twice as large as that used in 
previous calculations 0,E|. The number of radial functions per £ was 15 with the exception 
of £ = and £ = 1 where 17 and 16 functions respectively were used. The largest J was 13 
while L int was set to 4. The basis functions for each £ used a common exponent that had been 
energy optimized during some preliminary and smaller calculations. It must be emphasized 
that choosing a common A for both electron and positron states was not arbitrary but was 
a consequence of energy optimization process. The results of the CI calculation are listed in 
Table US 

The variation of pe and pr with J in Figure |H1 reflects the behavior seen in Figures ^ 
and El The values of p are smaller than the predicted asymptotic limits and seem to be 
approaching the correct value. Computational constraints mean that the dimensions of the 
radial basis, e.g. 15 e + and e~ LTOs per £, are smaller than those in the e + -H and e + Cu 
calculations. The PsH radial basis is further from convergence than the basis sets used for 
the equivalent calculations upon e + -H and e + Cu. 

Figure |H1 gives p values taken from the CI calculations of Saito Q. Saito used Natural 
Orbital techniques to reduce the dimension of the final diagonalization while using an orbital 
basis with J = 9 (this is about the same size as in j^|). Saito estimated the variational limit 
at each J and the Figure |H] curves were derived from this "full CI limit" calculation. Although 
the Saito curves have irregularities, they exhibit a p versus J variation similar to the early 
Bromley and Mitroy calculation 

Figure El shows the variation of (e)°° vs J. Once again, the 3-term series, Method III 
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FIG. 8: The exponents pe and pr as a function of J for 2 different CI calculations of PsH. Both 
the basis 1 and basis 2 with Lj nt = 4. The third data set was taken from the "full CI" calculations 
of Saito 




FIG. 9: The PsH binding energy (in units of Hartree) with respect to the Ps+H threshold of 
-0.750 Hartree as a function of J. The different curves use different algorithms to estimate the 
J — > oo correction as discussed in the text. The close to converged SVM energy |44| is shown as 
the horizontal line for comparison purposes. 

achieves its asymptotic value at the smallest value of J. Methods II and p av achieve their 
limiting values near J = 10. Methods I and GL again tend to underestimate (e) 00 while 
Method p overestimates (e)°°. 

The best CI estimate of the PsH energy is obtained by adding an Li nt correction of 
9.5 x 10~ 5 Hartree (the difference between the L int = 4 and L int = 9 energies |5j) to the 
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FIG. 10: The PsH annihilation rate ((r)°° in units of 10 9 s 1 ) as a function of J. The close to 
converged SVM annihilation rate is shown as the horizontal line. 

Method III e of 0.038923 Hartree. The resulting binding energy of e — 0.039018 Hartree is 
only 0.38 % smaller than the SVM e of 0.0391674 Hartree. 

Figure shows the annihilation rate versus J. The CI calculation does not converge 
to the SVM annihilation rate since a radial basis of 15 LTOs per I is simply too small. 
The i?p entries in Table HP reveal a 13.5% increase in Ar 9 between the basis 1 and basis 2 
calculations. The conclusions that can be drawn from Figure IIVI under such circumstances 
are somewhat limited. But Method III is again susceptible to the accuracy of the Ar J 
increments and again gives a final (T) 00 that is smaller than Method I. Methods II and 
p av again give final estimates of (T)°° that are close together. Method p gives the largest 
estimate of (r)°° while Methods I and GL give the smallest. 

Another problem with Method III arose from the usage of the Davidson method to 
perform the matrix diagonalization. This only gives values of (r) J that are stable to 6-7 
significant digits and this leads to the irregularities in the (T)°° evident in Figure ITU1 The 
problem of decreased precision when using the Davidson algorithm had been previously 
noted in CI calculations of helium and may be generic to iterative matrix solvers. 



IV. COMMENT ON THE SCALING OF THE ANNIHILATION RATE 

There is one class of system that has not been studied in the present work, namely the 
close to threshold scattering of positrons from atoms that can bind a positron. The behavior 
of the AZ j increments is complicated by a parametric dependence on the scattering length, 
A 0, 0, |48|. Once J is large enough to formally bind a positron, the magnitude of the 
scattering length decreases as J increases. Since Z eS oc A 2 , the decrease in A as J increases 
impinges on the increase in Z^ s that would otherwise occur. Indeed, one of the reasons why 
calculations on e + Cu were originally taken to J = 18 was to minimize the disruption that 
the scattering length had on Z e g pi. 

It was not worthwhile to try and analyze the behavior of the partial wave expansion for 
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FIG. 11: Plot of vs (J+i) 1 for the synthetic data as discussed in the text. The line shows the 
fit to eq. (|16j) while the • give the synthetic data points with the AZ^ S ~ ( J + I) -2 ' 6 dependence. 

the e + + Cu scattering system with its two major complications, the effect of the radial 
basis set and, secondly, the effect of the A versus J variation. Such an investigation is best 
delayed until substantially larger basis sets can be deployed. 

V. COMMENTARY ON THE WORK OF GRIBAKIN AND LUDLOW 



The problems caused by the slow convergence of Z e g can be exposed by a detailed anal- 
ysis of the recent Gribakin and Ludlow |14j calculation of the annihilation rate for e + -H 
scattering. This calculation used a variant of many-body perturbation theory (MBPT) to 
generate the initial set of (Z eS ) J data. They then did a fit to eq. ([To]) over the J = 7 — > 10 
interval and then plotted the results of that fit under the assumption that this demonstrated 
that their data obeyed eq. (fTE|). However, this procedure as applied by GL could hardly 
have been better suited to concealing deviations of the sequence of (Z e s) J values from the 
leading order (J + ±) -2 term (note, in this section the GL acronym refers to the results 
presented in I n effect, Figures 7 and 8 of GL are actually a demonstration of Taylor's 

theorem, namely that any continuous function will approximate a straight line if examined 
over a sufficiently small domain. 

To illustrate this point, consider a sequence of synthetic (Z cS ) J data generated by the 
prescription, (Z cS ) 7 = 2.5 and AZ^ S = 20 x (J + |) -2 ' 6 (such a data sequence could be 
generated exactly from a 3-term expansion of eq. ()14|) ). A fit to this sequence was made 
using eq. (|16|) and the results of that fit are shown in Figure ^2 Even though the data were 
generated according to p — 2.6, a visual inspection (and it should be noted that Figure [TT] 
is much higher resolution than Figures 7 and 8 of GL) suggests that the data are in good 
agreement with a, p = 2 power law decay whereas this is certainly not the case. 

The pz exponent derived from the recent GL calculation [14j is depicted in Figure [TJ 
Although GL assume pz = 2 when making the extrapolation, the actual exponent extracted 
from their sequence of Z e g values is p = 2.233 at J = 10. The present basis 2 calculation 
gives pz = 1.885 at J = 10. It is obvious that the GL calculation overestimates the rate 
at which AZ^ S increments are decreasing. This implies that the GL estimates of AZ^ S will 
be too small at higher J and this is the case. GL get AZ^ = 0.0501 while the CI-Kohn 
calculation gives AZ^ = 0.0604. Consequently, the GL calculation gives AZ^ 



;/>io = 0.48, 
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while the basis 2 CI-Kohn calculation gives AZ^ 10 ~ 0.62. The assertion by GL that their 
calculation has converged to the region in which the AZ^ S = B(J + §)~ 2 formula is valid is 
incorrect. 

Besides directly leading to less reliable extrapolation corrections, a related problem with 
the GL procedure is that it does not have the sensitivity to flag potential problems with 
the radial basis. A plot of pz versus J that crosses the pz = 2 line is a good indicator of 
some inadequacy in the basis. There is no indication that GL were aware that their AZ^ S 
were decreasing much too quickly as J increased, their statement that The use of a B-spline 
basis means that fast convergence is achieved with respect to the number of states with a 
particular angular momentum is difficult to reconcile with the present analysis. However, it 
should be noted that do indicate that they could improve the quality of their answers by 
"pushing harder" the numerics. 

The tendency for the GL calculation to overestimate the convergence of the annihilation 
rate increments probably does not arise from MBPT per se, rather it most likely comes from 
the underlying single electron basis. Besides the inherently slower convergence at higher 
J mentioned earlier, another possibility is due to the confinement of the basis to a box of 
radius 15 an- Confining the basis in this way will result in mean excitation energies (e.g. for 
predicting the multi-pole polarizab fifties) that will eventually increase as ~ £ 2 , where I is 
the orbital angular momentum, while for a real H atom the mean excitation energy for any I 
is less than 1.0 Hartree j4^|. Thus the occupancy of the higher J orbitals, which contribute 
significantly to Z e s will be inhibited, and successive AZ^ S will decrease too rapidly with 
increasing J. 

The relevance of these issues is best illustrated by a comparison with the exact value of 
Z e s which will be taken to be 3.407 at k — 0.4 a$ . This is 0.295 larger than the GL value of 
3.112. The underestimation of the higher partial wave contribution in the GL calculation, 
estimated at 0.14 = 0.62 — 0.48 is responsible for about 50% of the existing discrepancy. 

The discrepancies of GL with the best calculations are are not that severe for s-wave scat- 
tering since imposition of the p = 2 condition for J > 10 prevents the inherent deficiencies 
in their .B-spline basis from becoming too excessive. Also s-wave e + -H scattering system is 
certainly one of the easier positron annihilation calculations. However, the inadequacy of 
the GL methodology manifests itself more severely in other positron annihilation situations. 

The CI expansion converges quicker for electron-positron annihilations that take place 
at small distances from the nucleus than for annihilations that take place at large distances 
lEEEl- The presence of the centrifugal barrier for L > scattering leads to the electron- 



positron annihilations occurring further from the nucleus. Consequently the convergence 
problem is more serious for p-wave and d-wave scattering since a proportionally larger part 
of Z cS comes from the J — > oo correction [l^, EH- We have not repeated the earlier p- 
wave calculations jl^] with a larger radial basis, but a comparison with the CI-Kohn data 
at k — 0.4 <2q 1 indicates that the GL calculation again underestimates the impact of the 
high J orbitals. The CI-Kohn calculation reported in 2| gave AZ^ = 0.0480 while the 
GL calculation gave 0.0416. In addition, the CI-Kohn calculation gave p = 1.852 for J = 
10 while the GL calculation gave p = 2.152 (it is likely that an infinite basis CI-Kohn 
calculation would have p < 1.80 at J = 10). The GL calculation (which gave Z rme ff = 1.607) 
underestimates the L = 1 Z eS by about 0.18 at k = 0.4 an 1 (the T-matrix calculation gives 
1.786 BE} while van Reeth et.al. gave 1.794 QHEj). It is likely that at least 0.10 of 



the discrepancy will arise from orbitals with J > 10. 

One of the major results of the GL calculations was their demonstration that the en- 
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hancement factor is independent of energy. The enhancement factor can be defined as the 
factor that the annihilation rate, calculated as a simple product of the electron and positron 
densities needs to be increased in order to agree with the exact annihilation rate [13J, |48j, |50j ■ 
They (GL) based this conclusion solely on a forensic analysis of the Z e g annihilation rate 
matrix element. Figure 13 of GL reveals that the variation of the <i-wave enhancement fac- 
tor with energy is noticeably larger than either the s-wave or p-wave enhancement factor 
Since a larger fraction of the c?-wave Z e g comes from J > 10, the possibility exists 
that this stronger energy dependence is due to extrapolation issues as opposed to dynamical 
effects. Although GL seem unaware of the result, the slow variation of the enhancement 
factor with energy had been demonstrated in a model potential analysis [48|. Comparisons 
of model potential calculations with ab-initio variational and polarized orbital calculations 
had shown that a model potential calculation tuned to reproduce the energy dependence of 
the phase shifts also gave the energy dependence of Z eS [48]. The variation of the ci-wave 
enhancement factor was determined by tuning a model potential to the large basis phase 
shifts of [51] and then normalizing to a similar calculation of Z e $ The variation in the 
rf-wave enhancement factor over the energy range from k = to 0.5 1 was less than 4%. 
Although the model potential result is not conclusive, it does appear that the variation in 
the ci-wave Z e s is less than that indicated by the GL calculation. 

Another area where application of the GL methodology could lead to larger than an- 
ticipated errors is in the determination the angular correlation or the 7 energy spectrum 
[o^ . These two properties depend on the relative momentum of the annihilating electron- 
positron pair [45J, |53j . It is known from investigations of momentum space wave functions 
that the low momentum part of the wave function largely arise from the large r part of the 
wave function while the high momentum properties come from the small r part of the wave 
function jH^. Under such circumstances, application of the GL method could easily result 
in errors to the J —>■ 00 corrections that depend systematically on the 7-energy or recoil 
momentum. 



VI. SUMMARY AND CONCLUSIONS 

Single center methods represent a superficially attractive method to study mixed electron 
positron systems since existing computer codes can be adapted without too much effort. The 
penalty associated with this approach is the slow convergence of the binding energy and, 
more noticeably, the annihilation rate with respect to the partial wave expansion of the single 
particle basis. The results presented here are generally consistent with the asymptotic limits 
derived from second-order perturbation theory by Gribakin and Ludlow |28j |. The actual 
calculations at finite J generally give increments to the energy (phase shift) and annihilation 
rate that decrease slightly slower than the GL limits, but the overall trends are compatible 
with the GL limits. 

The tendency for the convergence with respect to the radial basis size to slow down as 
J increases does have implications for the design of any CI type calculation. Some sort of 
extrapolation in J is necessary in order to determine the energy and more particularly the 
annihilation rate. But there is no point in making J bigger if this is done at the expense of 
the radial basis set. One simply ends up with increments to the energy or annihilation rate 
which are systematically too small at higher J. This problem does not seem to be restricted 
to the LTO basis used in the present work. Convergence problems at high J are also present 
for the Gribakin and Ludlow calculations which used a S-spline basis [l4| and the Saito 
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calculations which used a natural orbital basis 22|. It is amusin g to note that the one 

of the first manifestations of this problem occurred over 40 years ago |29l . 55l . 

The best methods for estimating the J — > oo corrections depends on the quality of 
the underlying calculation. For a low precision calculation, Method I would seem to be 
appropriate. A low precision calculation can probably be regarded as one with p E or p r 
exceeding 4 or 2 respectively when the J — > oo correction is evaluated (assuming that p 
approaches its limiting value from below). Methods II or III would seem to be the preferred 
options for a high precision calculation. As a general principle, inclusion of the second term 
in the asymptotic series leads to improved (X)°° predictions when compared with asymptotic 
series with the single term series. Method III is more susceptible to imperfections in the 
radial basis and should not be applied to the calculation of the annihilation rate unless a 
very large radial basis set is employed. Irrespective of how the J — > oo corrections are 
evaluated, it is essential that the exponents p relating the changes in the expectation values 
be examined test of the quality of the radial basis. 

The overall situation regarding the use of single center methods to compute positron- 
atom phase shifts or energies is that calculation to the sub 1% accuracy level is achievable 
for those systems that have a parent atom ionization potential greater than 0.250 Hartree. 
The 0((J+ §)~ 4 ) convergence means a J of 10 or slightly larger will generally suffice as long 
as the method used to perform the J — > oo correction is more sophisticated than those used 
previously. However, the situation with respect to the annihilation rate is much grimmer 
and it is not possible to guarantee 1% accuracy for even the simple e + -H system. Here the 
0((J + |) -2 ) convergence means the J —>■ oo correction is larger, and moreover the slow 
convergence with respect to the radial basis is further complicated by the fact that it is 
slower at high J than low J. In this case, it appears that "God is on the side of the big 
basis set" 15611 . 
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TABLE III: Results of CI calculations on e + Cu versus J for the energy and annihilation rate for 
a series of J. The binding energy of the positron to neutral Cu is denoted e (the energy of Cu 
with respect to the core was -0.28394227 Hartree Q). The spin-averaged 27 annihilation rates are 
given for the core (T c ) and valence (T v ) electrons. The results in the row 00 include the J — > 00 
correction evaluated at J = 18 using the various methods described in the text. 



J 


£ 


r c 


r„ 




Hartree 


10 9 sec' 1 


10 9 sec" 1 





-0.00112467 


0.000289 


0.000132 


1 


-0.00080292 


0.001443 


0.001692 


2 


-0.00037356 


0.004818 


0.009728 


3 


0.00031179 


0.011213 


0.033656 


4 


0.00111995 


0.017605 


0.070360 


5 


0.00187958 


0.022223 


0.110121 


6 


0.00251736 


0.025271 


0.147793 


7 


0.00302852 


0.027272 


0.181727 


8 


0.00343136 


0.028611 


0.211657 


9 


0.00374774 


0.029526 


0.237830 


10 


0.00399692 


0.030168 


0.260656 


11 


0.00419429 


0.030627 


0.280571 


12 


0.00435174 


0.030963 


0.297984 


13 


0.00447832 


0.031212 


0.313256 


14 


0.00458086 


0.031402 


0.326694 


15 


0.00466456 


0.031547 


0.338564 


16 


0.00473338 


0.031660 


0.349087 


17 


0.00479037 


0.031749 


0.358454 


18 


0.00483788 


0.031821 


0.366821 




J — > 00 






P 


3.2751 


4.0511 


2.0295 


Met hot I p 


U.UUozUlz 


U.UoZZiy 


U.oloU ( 


Method p av 


0.0051481 




0.51526 


Method GL 


0.0050998 




0.51325 


Method I 


0.0051080 




0.51751 


Method II 


0.0051584 




0.51529 


Method III 


0.0051603 




0.49364 


Method S 


0.0051394 




0.46398 




Earlier calculations 






FCSVM [8, 43] 


0.005597 


0.0339 


0.544 


CI, J = 18 [8J 


0.004786 


0.03173 


0.35499 


CI, J -» 00 [8J 


0.005117 


0.0321 


0.4744 
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TABLE IV: Results of CI calculations on PsH for orbital bases with L- Uit = 4 and for a series of 
J. The total number of electron and positron orbitals are denoted by iV or b, the total number of 
configurations is given by Nqi, and the LTO exponent for i = J is listed in the A column. The 
3-body energy of the PsH system in Hartree is denoted by E^PsH) and the T is given in 10 9 sec . 
The J — > oo extrapolations were carried out at J = 13. 



J 


A 


N orh 


Nci 


£(PsH) 


r 





2.10 


17 


2601 


-0.69133618 


0.374196 


1 


2.26 


33 


9265 


-0.74705969 


0.782256 


2 


2.36 


48 


22810 


-0.76620031 


1.080456 


3 


2.46 


63 


44650 


-0.77514128 


1.292538 


4 


2.52 


78 


78640 


-0.77995286 


1.448216 


5 


2.72 


93 


120265 


-0.78274494 


1.566206 


6 


2.93 


108 


165265 


-0.78449597 


1.658344 


7 


3.13 


123 


213415 


-0.78565220 


1.732023 


8 


3.34 


138 


263365 


-0.78644538 


1.792061 


9 


3.56 


153 


314890 


-0.78700639 


1.841756 


10 


3.75 


168 


366415 


-0.78741330 


1.883366 


11 


3.95 


183 


417940 


-0.78771485 


1.918595 


12 


4.15 


198 


469465 


-0.78794247 


1.948689 


13 


4.35 


213 


520990 


-0.78811707 


1.974632 






J - 


-> oo limits 






P 








3.44547 


1.92375 


Method p 








-0.7889974 


2.3352 


Method p av 








-0.7888995 


Z.oZoZ 


Method GL 








-0.7887894 


2.2988 


Method I 








-0.7888198 


2.3121 


Method II 








-0.7889218 


2.3225 


Method III 








-0.7889231 


2.2915 


Method S 








-0.7888843 


2.2203 






Other calculations and earlier CI calculations 




SVM [44] 








-0.78919674 


2.4712 


J = 9, Li nt = 


4 [51 


90/91 


63492 


-0.7866818 


1.7903 


J = 9, -C/int = 


9 [5J 


90/91 


95324 


-0.7867761 


1.7913 


J = 9, [9J 




oc 


oc 


-0.786949 


1.8230 
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